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Abstract. We determine the three loop anomalous dimensions of the quark, centre and off- 
diagonal gluons, centre and off-diagonal ghosts and the gauge hxing parameters in the maximal 
abelian gauge for an arbitrary colour group in the MS renormalization scheme at three loops. 
We show that the three loop MS /5-function emerges from the renormalization of the centre gluon 
and also deduce the anomalous dimension of the BRST invariant dimension two mass operator. 
Moreover, we demonstrate that in the limit that the dimension of the centre of the group tends 
to zero, the anomalous dimensions of the quarks, off-diagonal gluons and off-diagonal ghosts 
tend to those of the quarks, gluons and ghosts of the Curci-Ferrari gauge respectively. 
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1 Introduction. 



The multiloop renormalization of quantum chromo dynamics (QCD), the quantum field theory 
underlying the strong interactions, has now been successfully determined at four loops in the 
MS scheme, El El EJ El E3 Hj • Indeed the one loop ^-function, jQj, establishes the important 
property of asymptotic freedom. Further, with the need for more accurate theoretical results 
such as the precise way in which the coupling constant runs, higher loop corrections proved 
necessary. Subsequently, the scheme independent two loop result was computed in [2] prior to 
the three loop calculation of 4 . Given the large increase in the number of Feynman diagrams 
with loop order and the parallel problem of devising an algorithm to extract the divergence 
structure of difficult four loop master integrals, it was several years before the four loop (5- 
function appeared, 6 ]• Indeed given the complexity of such a calculation, it was only technically 
possible with the intense use of the symbolic manipulation programme Form, 8 . Though the 
three loop result of [Ij also used computer technology and the Mincer algorithm, There was 
an underlying thread to all these computations which lay in a judicious choice of gauge in which 
to perform the calculation. Although the /3-function is gauge independent, choosing a general 
covariant gauge, say, to carry out the calculations could have resulted in a large amount of extra 
unnecessary computation. This was avoided by considering the Feynman gauge where the gluon 
propagator reduces to one term proportional to a scalar field propagator. Only after the original 
Feynman gauge calculations were performed were computations with gluon propagators in the 
full covariant gauge subsequently carried out, 00 [7j. These were necessary for other problems 
aside from justifying the full gauge parameter independence of the /3-function. 

For instance, the anomalous dimensions of the fields as functions of the covariant gauge 
parameter, a, were required for a variety of composite operator renormalizations such as those 
central to deep inelastic scattering. (See, for example, ^^.) Also, it has recently been 
established that there is an interesting relation, |121 I13| . in respect of the dimension two BRST 
invariant operator which could play the role of a gluon mass. In |12| 113] it was demonstrated 
that in the Landau gauge its renormalization is not independent, being related to the gluon 
and ghost anomalous dimensions. This was observed by an explicit three loop computation 
in the MS scheme, |12| . More recently, the explicit renormalization has been determined at 
four loops through the provision of the Landau gauge gluon and ghost anomalous dimensions 
at that order, ^Ij. Significantly, similar identities for the analogous operator exist in other 
gauges such as the maximal abelian gauge (MAG), ^3E3> and in space-time dimensions other 
than four, [Uj. Since these dimension two operators have been the subject of intense analytic 
investigation in various gauges in recent years, see, for instance |T§1 El HOI [2D HH HHI 123 HH 
l26| |23 EHl EH1 EH EH an d references therein, due to their condensation in a non-trivial 
vacuum, there is a clear need to renormalize QCD in this gauge. In particular the explicit 
values of all the anomalous dimensions are required as the first step in the extension of the 
local composite operator (LCO) method for QCD, ^H]; to the MAG in various colour groups. 
This would thus open up the possibility of extending the effective potential calculations in the 
Landau gauge, ^Bl EHl- to situations beyond the few one loop SU(2) MAG studies already 
considered, |19l 1221 1261 1321 133j . This is the main aim of this article where we will perform the 
full MS renormalization of QCD in the MAG for an arbitrary colour group to determine the 
explicit values of the anomalous dimensions with the renormalizability of the gauge having been 
discussed in (221 EU EH EH EH EE! ■ Though given the nature of the MAG construction where the 
colour group is split into its centre and off-diagonal sectors, we will make several assumptions 
about the group structure which we have checked are at least valid in SU(2) and SU(3). It is 
important to note that the only previous explicit renormalization of QCD in the MAG was at 
one loop and for the specific group SU(2), [23 E21 E3 EH] • 
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In referring to the MAG it is important to note at the outset that we are in fact considering 
the more general modified MAG as discussed in j^2] for SU(2). The reason for this is that 
the true MAG is defined in a similar fashion to the Landau gauge. However, by minimizing the 
square of the gauge potential over only the off-diagonal sector of the the colour group, as opposed 
to the full group in the usual covariant gauge situation, it transpires that the renormalization 
of the subsequent gauge fixed Lagrangian is singular. Therefore, analogously to the generalized 
Landau gauge or covariant gauge, a covariant gauge parameter, a, is introduced which is not 
to be confused with the parameter of the covariant gauges. With this non-zero a one has the 
modified MAG and as we will show, it is the renormalization of a itself which becomes singular 
as a — > 0. However, all the remaining renormalization group functions are finite as a — > whence 
one obtains the true MAG anomalous dimensions. Moreover, as has been observed before, |19j . 
the structure of the MAG renormalization has connections not only with the Landau gauge but 
also with the related non-linear covariant gauge known as the Curci-Ferrari gauge introduced in 
[HH] . It will turn out that such connections will also prove useful for justifying our final three 
loop MS anomalous dimensions. 

Another motivation for considering the MAG rests in one of the original reasons why it was 
introduced. One possibility for the mechanism of confinement is the condensation of abelian 
monopoles which clearly originate from the centre of the colour group, |4L)l 1411 W2\ . In any cal- 
culations which seek to focus on this supposition, it makes sense to consider a gauge where the 
centre and off-diagonal fields are separately identified in the gauge fixing. Therefore, by estab- 
lishing the renormalization structure at three loops in this gauge, one would expect the results 
will be useful, say, in any continuum matching one might have to do in lattice computations. 
On a final note we draw attention to another gauge in which QCD is renormalized and that 
is the background field gauge where the gauge field is split into a classical and quantum part, 
[131 El G3 HE] • The latter is regarded as the totally internal quantum fluctuation. In addition 
to the other three loop results referred to earlier, QCD has also been renormalized to the same 
order in this gauge, |451 1461 H7] . The main advantage of the background field gauge is the fact 
that the f3- function emerges from the renormalization of the gluon field. In other words one 
needs only to consider a 2-point function rather than a 3-point function which considerably sim- 
plifies any explicit computation. Interestingly, the MAG, where the gluon field is split, but with 
respect to the colour property, has an analogous simplification which is that the centre gluon 
anomalous dimension is also equivalent to the /^-function, jS2]. This feature will be exploited 
here to reduce the number of Feynman diagrams we have to consider to perform the full three 
loop renormalization. 

The paper is organized as follows. In section 2 we review how the MAG Lagrangian itself 
is constructed prior to summarizing the group theory results which were required for the three 
loop renormalization. This is a non-trivial exercise since the colour indices have to be identified 
either as originating in the centre of the Lie group or in the off-diagonal sector. The details 
of the full three loop renormalization are discussed in section 3 where the structure of the ac- 
tual renormalization established with the algebraic renormalization formalism, |32j . is reviewed. 
This section also contains the main results of the computation which is the determination of 
the explicit values of all renormalization group functions for the MAG. Finally, section 4 con- 
tains concluding remarks and the appendix contains the non-trivial Feynman rules used in the 
calculation. 
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2 Maximal abelian gauge. 



We begin by recalling the essential features of the maximal abelian gauge fixing which depends 
on the parameter a. First, we note that the colour group generators are T A where 1 < A < Na 
and Na is the dimension of the adjoint representation. Thus the group valued gauge field A M 
can be decomposed as 

A, = A A T A . (2.1) 

In considering the MAG the group generators are split into two sets. Those corresponding to 
the generators of the centre of the group, which themselves form a group, and the remaining 
set. For notational purposes we will use the indices i, j, k and I to denote centre elements and 
a, b, c and d to denote off-diagonal elements. Thus A M can alternatively be decomposed as 

A? = A a ^T a + A\T (2.2) 

where we introduce the dimension of the centre by noting that 1 < i < N^ and allowing the 
off-diagonal indices to range over 1 < a < N£. Clearly 

N£ + NX = N A (2.3) 

and, for instance, in the unitary groups SU(N C ) we have N$ = N C -1 and iVf = N c (N c -l). With 
this notation the QCD Lagrangian in general is, with the gauge fixing part Lgf to be specified, 

L = — \g%G a ^ + i^W + % (2.4) 

where G A = d^A^ — d v A A + gf ABC A^A^, is the covariant derivative, there are Nf flavours 
of quarks, Np is the dimension of the fundamental representation and g is the coupling constant. 
For the MAG the indices A are split into the two sectors giving 

L = — \G%G a ^ - ~G^G^ + W + L gf (2.5) 

where now Lgf is interpreted as the MAG gauge fixing term. This is constructed, see, for 
example, |22l I32j . in the standard way by the BRST variation of a specific operator. In the 
usual covariant gauge fixing one uses 



Lgf = 55 



lA A A AfJ " + \ac A c A 



(2.6) 



where 5 and 5 are the BRST and anti-BRST variations respectively, c A is the ghost field and 
c A is the anti-ghost field. In the MAG the gauge fixing term is chosen in a similar way. The 
off-diagonal sector is chosen as in the covariant gauge case but the diagonal sector is restricted 
to being in the Landau gauge to fully fix the gauge overall. It is not instructive to repeat all the 
additional technical details of the gauge fixing which have been discussed previously in [221 02] • 
Therefore, for the MAG we take, [221132] . 



Lgf = 55 



l A a A a„ + i a gi c a + l^^M + (!_£)$ ^AJ, (2.7) 



where the last term is included to ensure one can interpolate the results between the MAG 
and the Landau gauge according to how one chooses the additional parameter For instance, 
the Landau gauge corresponds to a = and £ = 1 and the (modified) MAG is a / but 
£ = 0. This particular gauge fixing was introduced in [22] and we have chosen to work with 
this version for various reasons. First, this Lagrangian has been examined from the algebraic 
renormalization point of view and the Slavnov- Taylor identities have been established. Second, 
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and more crucially for the current article, in a computation of the magnitude of the three loop 
MAG renormalization it is important to recognise that calculating with an arbitrary a and 
£ allows us to check the correctness of, say, our programming and resultant renormalization 
constants. In particular the Landau gauge three loop anomalous dimensions ought to correctly 
emerge from the computation prior to specifying the MAG values of the parameter £. This is 
actually a non-trivial point since we have to perform the group theory manipulations for the split 
group and not the full group as one would do in an ordinary covariant gauge fixed Lagrangian 
where the Casimir structure resulting from group identities is already well established. 

With the MAG gauge fixing, (|2.6|) . it is elementary to perform the BRST and anti-BRST 
variations, which are given by 



5c a 



(<V a + 9f ajc Aic c + gf abc Ay + gf abk Ay 



g f abk c b c k + |/ a ~f r 



eabcb c 



dy + gf ibc Alc c 



5(f 
5c l 



If 
1 



5b a = , 



ribc b c 

9 7 c c 



bb l 



(2. 



and 



5c a 
5c a 

6 A 



(d,c a + gr c A^c c + gf abc A b ^ + gf abk Ay 



b a + g fabc c b^ + g fabk c b-k + g jabkJ>Js 



C C 



gf abk-b-k + LgfabcJ^ 



I' 

5b { 



8b a 
8<? 



gf abc b b^ _ g fa b%b -K + g fa0^ b K ^ 

1 



pabkrb-k 



cabk-hik 



- gf ibc b b (? 

respectively, to obtain the gauge fixed MAG Lagrangian 



V + gf lbc c b s 



Sc l 



-gf lbc c b c c 
2 



(2.9) 



-gb a 

+ g 



l 



(1 - ()f abk A b A k v - -af abc c b c c - af abk c^c 



(1 - ()f abk A a c k d^c b - (f abk A a d> M c b c k - f abc A a d tl c b ? 



\o.^t> „k 



^abk A a bg^-k 



- f abk d'"A k ^c b 

(1 - Or d cbd AlA b ^c d + (1 - C)fo dcj AlA jl *(?c d + (1 - ois lcj K Aj ^ cl 



+ (l-{)fV d %A^c?c d -jf t 



a -abcd^L^) c c c d 



a^bcd^a-b c c c d 



jf* bd (?c b c c c l 



(2.10) 



where we have introduced the term \ab l b l to fix the residual gauge freedom in the full gauge 
field, 22- We have also simplified the notation by defining 

iAB riCD rABCD _ reAB reCD 



ABCD 
d 



J o 



2.11 



for the 4-point vertices. For perturbative computations the auxiliary fields b a and b l are elimi- 
nated by their equations of motion 



b' 



1 

a 
1 

a 



d»Al - (1 - C,)gf ahk AlA k » - -agf abc ^c c - agf abk c b c k 

d^Ai 



(2.12) 
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to obtain the MAG Lagrangian in the form we will renormalize it, |32| . 



+ 9 



[l - Q)f abk Alc k d»c b - Qf abk Ald»c b c k - f abc A a Jd^c c - Cf abk A a Jd^c k 



(i-C) 



f abk d^A a tl A b v A kv - f abk d^A a llC b c k - -f abc d^A^c c 
- (2 - C)f abk A k c a d >1 c b - f^d^A^c 

(1 - 0f2 chd AlA b ^(fc d - ^~^ 2 jf bl AlA bil A k v A lv + (1 - Qfo dcj KA j ^c c c d 



+ 9 Z 



2a 



(i-C) 



fWed A a A j^ c d + ^ _ Qf^cl A a A j^ c l + ^ _ Qf^A%A jl *Fc l 



(1 - 0fo jdi AlA j v<?c d - jff cd (f(?c c c d 



-ff cd ^^c c c d + -f% cbd ^c c c d 



-jf bc c a ^c c c l + -ff bl ^c b c k c l 



(2.13) 



where it is understood that the parameter a, which is distinct from a, is set to zero after our 
renormalization. Having constructed the full MAG Lagrangian as a function of the parameters 
a, C, and a we note that the full set of non-zero Feynman rules generated from (|2.13jl are given 
in appendix A. 

Since we will be performing our calculation for an arbitrary colour group but with the group 
algebra split into centre and off-diagonal sectors, we close this section by discussing the main 
properties of the Lie algebra which were required. To construct the necessary lemmas we recall 
that the Lie algebra and basic Casimirs for the full group as well as the Jacobi identity are 



rpA rpB 



ACD f BCD 



ABCrpC 



if 

C A 5 AB 
f AB V 



rpArpA = 

ABE rCDE . f BCE f ADE 



, Tr(r A r B ) 

_j_ j-CAE jBDE 



T F 5 AB 



(2.14) 



From the Lie algebra we have that = and f l ^ a = which enshrines the centre property 
in the algebraic manipulations. So the second equation of Q2.14JI gives the relations 



C A 5 ab 
C A S ij 



jacd jbcd _|_ 2 j ac 3 j^bcj 
jicd jjcd 



(2.15) 



To proceed we make the assumption that f acd f bcd is proportional to 5 ab which is certainly true 
for SU(2) and we have checked it is also valid in SU(3). In the group theory discussion which 
follows, it is important to bear in mind that for groups where this simplifying feature is not 
present then one would have to proceed with f acd f bcd being proportional to a symmetric rank 
two tensor. Taking contractions of Q2.15|) leads to 



jiab jiab 



n£c a , f abc f abc 



NX - 2N% 



C A 



(2.16) 



where is the dimension of the centre and is the dimension of the complement of the 
centre. Hence, 

jicd rjcd _ (J A ff3 



jacj jbcj 



N d 

A r~i tab facd fbcd 

nx Ca6 ' f f 



[NX-2Nj] b 



(2.17) 
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With these elementary results we can use the Jacobi identity to establish several useful relations 
which were used extensively throughout the computation 



japq jbpr jcqr 
jipq jbpr jcqr 



MJ (j^jabc japq jbpi jcqi 



2NX 



C A f 



(the 



[NX - 2N$\ 



2NX 



C A f ibc fiPQ fbpj -P 



<:<!.) 



Nl 
2NX 

A f~\ fibc 
NX ' 



For the group generators, we have 



Tr(T a T fe ) = T F 5 ab , Tr (r a T') = , Tr (T'TA = T F 5 lj 



(2.18) 



(2.19) 



and we make the assumption that T l T l is proportional to the unit matrix which is certainly true 
for the groups SU(N). Then, from (|2.14|) we have 



rplrpl 



Tf 

N F ' 



mi 



after contracting Tr (T l T J ). Hence, 

rparpd 



I . 



It therefore follows from the Lie algebra itself that 

C A 



C F - 

—NX 

N F A 



Tf d C A Ni 
Nf Na + ~2NX 



rj-il rj-iCL rpl 



Tf 
N F ' 



Na 



C A Nj 
2NX 



Ca 
2 



rpjrpirpj _ Ardrpi 

N F A 



(2.20) 
(2.21) 

rj-iCL 

2.22) 



As a consistency check on these results adding the first pair together recovers the usual result 



C F 



C a 



(2.23) 



for a free off-diagonal index. Summing the final pair is also consistent with this result after use 
of the relation 



C F N F 



N° A + Na 



T f 



(2.24) 



which follows from taking the trace of T A T A . Next, given the Lie algebra it is straightforward 
to construct the useful lemmas 



jabcrpbrpc 

jibcrpbjnC 



i[NX-2Nj] 
2NX 

\CaT. 



C A T a , jabjrpbrpj 



iNi 

2NX 



(2.25) 



Whilst these results proved to be the workhorse for the full three loop computation as well, it 
turned out that at three loops it was quicker to include additional lemmas to speed up the group 
theory computation of our Form programmes. These were derived from several applications of 
the Jacobi identities and are 

japq jbrs jqms jemt jprt q 

japq jbjs jqms jemt jpjt q 
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japq jbrs jqjs j c jt j'Prt q 

fapq fbrj fqmj fcmk fprk _ A A fabc 

J J J J J 4iV°^ 
fapj fbrs fjms fimt fprt N aI N A ~ 2N a\ C A f abi 

j j J j j 4i\T°^ 

N d2 C 2 

japj jbks jjms jimt jpkt A A jabi 

N A 

japq jbrs jqms jimt jprt q 

N d2 r 2 

japj jbks jjms jcmt jpkt A A jabc (2 26) 

where we note that the indices m, p, q, r, s and t are also regarded as off-diagonal. Given 
the structure of the three loop Green's functions these relations were sufficient for handling 
the group theory associated with the gluon 2-point functions. In that case for any Feynman 
diagram one has at most six structure functions contracted together with two free external group 
indices. However, in the renormalization of the A^c 1 ^ vertex at most seven structure functions 
are contracted together with three free group indices. For this case the Green's function was 
multiplied by an additional structure function to leave a scalar group string to be simplified. 
The route to achieving this, aside from applying the rules discussed so far, was to follow the 
approach used for the structure constants of the full group in that in that case they correspond 
to the group generators in the adjoint representation. In other words one replaces the structure 
constants by 

ftJu) BC = - *f ABC (2-27) 

and then applies the usual Lie algebra properties to T A ^ with the proviso that one evaluates 
identities in the adjoint representation. For instance, the result 

t b t c t a t b t c = [c f -Ca){C f -\C a )T a (2.28) 

implies 

r 

1 adj 1 adj 1 adj 1 adj 1 adj 

For the MAG calculation one can also use this strategy provided one appreciates that the struc- 
ture constants with two or more centre indices are identically zero which means that the non- 
trivial structure constants have at least two off-diagonal indices. These are therefore regarded 
as the matrix indices which leaves the third index as the free generator index and this can either 
be centre or off-diagonal. More significantly one must be careful in regarding these objects as 
nothing more than matrices and not as representations of the group generators since the matrix 
indices are only elements of the off-diagonal sector which is not closed in the group sense. In 
light of this to differentiate from the adjoint representation of the full group we therefore choose 
to define the analogous object S A by 



rpB rpC rpA rjiB rpC f\ ( O on\ 

i adi i adi i adi i adi i adi ~ U ■ (Z.ZV) 



(S a ) bc = f abc , (S l ) bc = f bc (2.30) 

where the two matrix indices will always be off-diagonal. Subsequently, given all the previous 
lemmas it only remains to resolve objects of the form 



and 



tr (S A S B S C S D ) tr (s A S B S c S D ) (2.31) 
tr (S A S B S C S D S A S B S C S D ) (2.32) 
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where the trace is the usual matrix trace and we use tr in contradistinction to the TV of the 
full group. Such structures are known to occur in higher loop calculations in QCD when S A is 
formally replaced by T A , but not until four loops where they arise only in terms involving 
the simple pole in e where d = 4 — 2e. Therefore, to ensure renormalizability they either have to 
vanish or cancel since at three loops they can potentially occur in the double and triple poles in 
e as well as the simple one. As these two structures emerge with the summed indices in various 
combinations of centre and off-diagonal indices, it is appropriate to relate them to a common 
term via the relations 



tr 



tr 



^ S i 3 j s k s 

(s i s j s k s d 



tr (S i S j S k S l 

(s i s j s k s d 



tr 



tr 



tr 



(s i s j s k s i 



tr (S i S j S k S d 



+ 



6N$ + NX 



4Nf 



tr 



Am 



o2 



A 



tr (s i S i S c S d ) tr (s i S i S c S d 



8NX 3 



- tr (s i s b s°s d ) tr (s i S b S c S l 

+ tr (s i S b S°S d S i S b S c S d ) + tr 



Id* 



5NX -4N£N£ + N£ 



o'2 



8N£ 3 



and 



tr(s i S j S c S d S i S j S c S d ^) = 



S c S c S*j 

(2.33) 

(2.34) 



which are readily established by use of the Lie algebra and the Jacobi identity. In the ac- 
tual renormalization of the A^<?c 6 vertex this leaves the two as yet unevaluated structures as 

tr (s i S b S c S d> j tr (VS^S^) and tr (s i S b S c S d S i S b S c S d ) . It turns out that when the pole parts 
of all the Feynman diagrams for this renormalization are added up then the coefficients of these 
structures is finite. 



3 Renormalization. 

Having derived the MAG Lagrangian we now turn to the details of its renormalization. First, in 
renormalizing a renormalizable quantum field theory one ordinarily introduces renormalization 
constants for all the fields and parameters in the Lagrangian. For field theories possessing 
symmetries such as a gauge symmetry these renormalization constants are not necessarily all 
independent. The underlying symmetry can constain several or more to be related. To determine 
such relations, one can apply techniques such as algebraic renormalization, |3E], which ensures 
the Lagrangian is stable under quantum corrections. In [^1] this approach has been applied to 
(|2.13|) and several interesting relations emerge. For instance, it turns out that the anomalous 
dimension of the centre gluons is proportional to the QCD /3-function. This is a useful result 
since for this gauge fixed version of QCD it means that one does not have to renormalize a 
3-point vertex to determine the known three loop /3-function of jlj. Instead one needs only to 
consider the centre gluon 2-point function. From a practical computational point of view this is 
a significant observation which we exploit later. Moreover, a similar property is also present in 
the background field gauge where the anomalous dimension of the background gluon is simply 
related to the coupling constant renormalization, |45[ I46j . Although we are a priori aware of the 
relation of the centre gluon renormalization to that of the coupling constant renormalization in 
the MAG, in defining our renormalization constants we choose at the outset to leave this result 
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to emerge in the computation rather than put restrictions on the initial setup. Therefore, we 
define the renormalization constants as 



A a » = JZ~ A A a » , 4" = y/Z^A*" , c a Q = ^Z c c a , cf Q = y%c? , 




(3.1) 

where \x is the renormalization scale introduced to ensure the coupling constant is dimensionless 
in d dimensions, the subscript denotes the bare quantity and the subscript i in the field 
subscripts of the renormalization constants is included to indicate that they correspond to centre 
objects and there is clearly no summation over repeated indices in this instance. In writing 
down (|3.1j) from 32| we have chosen, by contrast to A 1 , to encode the structure of the centre 
ghost renormalization. In particular the anti-centre ghost and centre ghost renormalizations 
are, contrary to the usual covariant gauge ghost renormalization, inverses of each other and not 
equal. This property emerges from the algebraic renormalization analysis, [H2]- Prom a practical 
point of view this means that the centre ghost 2-point function cannot be used to determine 
Z c i since it would be finite, |H2]- Instead to find Z c i one has to renormalize the 3-point A^c 1 ^ 
vertex once the coupling constant and off-diagonal gluon and off-diagonal ghost wave function 
renormalizations have been determined at that particular loop order. Therefore, in the MAG 
one has still at least one 3-point function renormalization to perform. 

However, the benefit in determining Z c i rests in the fact that the dimension two BRST 
invariant operator 

O = \AlA a » + ac a c a (3.2) 

possesses an interesting renormalization structure, ^3 EH- It transpires that its anomalous 
dimension is not independent but satisfies 

lo(a) = - ^ + 7c* («) (3-3) 



where 
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16tt 2 

and for completeness its associated renormalization constant is defined as 



(3.4) 



= Z o . (3.5) 

Therefore, it will be straightforward to deduce 70(0) from explicit knowledge of 7 c »(a). This is 
one of the key results required for a two loop extension of the LCO method to the condensation 
of O in the MAG and will be one of the main results of the article. That such a relation is 
present in the MAG is not specific to this gauge. A similar relation exists in the Landau gauge, 
|12| I13j. for the analogous operator where the indices range over the full colour group. We have 
also introduced a renormalization constant in (|3.1|) for the interpolating parameter £. However, 
since we are only interested in the renormalization of the MAG itself which corresponds to the 
fixed point value of £ = it turns out that for the MAG renormalization the explicit form of 
Zc is not required since it will always be multiplied by zero and there are no singularities in £ 
in the Feynman rules. 

We now turn to the technical details of the renormalization of (|2.13j) at three loops in the 
MS scheme. First, the renormalization group functions we will determine are deduced from the 
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Green's function 


One loop 


Two loop 


Three loop 


Total 


Aa Ab 


6 


131 


6590 


6727 


A 1 A> 


3 


54 


2527 


2584 


c a c b 


3 


81 


4006 


4090 


1p if) 


2 


27 


979 


1008 


A*c i c b 


5 


287 


22621 


22913 


Total 


19 


580 


36723 


37322 



Table 1. Number of Feynman diagrams for each Green's function for the MAG renormalization. 



explicit respective renormalization constants themselves for the MAG, via 

7a(<0 
la (a) 
7 A* (a) 
la- (a) 
lc (a) 

7c* ( a ) 
7v( a ) 
70 (a) 

where, similar to the Curci- Ferrari gauge, we have not assumed that Z a = 1. Though we 
have set a = and £ = 0. Since these are the renormalization group functions we require, we 
will therefore renormalize the centre and off-diagonal gluon 2-point function, the off-diagonal 
ghost 2-point function, the quark 2-point function and the A^cfc 6 3-point function all at three 
loops. For these Green's functions the Feynman diagrams were generated with the Qgraf 
package, [IH], and the specific number of diagrams at each loop order and Green's function are 
summarized in Table 1. By contrast, to indicate the magnitude of the MAG renormalization 
using the Feynman rules of the appendix, we have provided a similar diagram count for the 
Curci- Ferrari gauge three loop renormalization of ^2] in Table 2. To proceed we convert the 
Qgraf output format to the electronic notation used by the Mincer algorithm, j^j, as written 
in the symbolic manipulation package Form, (HUSH) in terms of diagram topology and internal 
momentum routing. The Mincer algorithm is then applied to all 37322 Feynman diagrams 
required for the full renormalization. Though it ought to be noted that given that the main 
group theory is carried out prior to determining the divergence structure of a diagram the value 
of a graph could be zero purely from group considerations. For instance, when one has a one 
loop gluonic self-energy subgraph anywhere where one external leg of the subgraph is in the 
centre of the group and the other in the off-diagonal sector, then that graph is trivially zero 
since f lcd f bcd = o. To appreciate the benefit of the centre gluon anomalous dimension relation 
to the /3-function, if such a relation did not exist then one would have to compute a 3-point 
function in addition to the ones listed in Table 1. The easiest one from a computer algebra point 
of view is the quark gluon vertex. For an off-diagonal gluon the figures for the corresponding 



d d 
(3(a)— \nZ A + ai a (a) — In Z A 



d 

(3(a)— \nZ a - ia{o) 



1 - a— lnZ a 
oa 



(3(a 
(5{a 

f3(a 
(3(a 



J^lnZ^ + ai a (a) ^ In Z Ai 

d d 

— In Z a r + ai a (a) — In Z al - 7 A i (a) 
oa aa 

— lnZ c + ai a (a) — In Z c 

tt- In Z c i + ai a (a) In Z c i 
oa oa 

— InZ^p + ai a (a)—\nZ^ 
J^lnZo + aia(a)-^lnZ 



(3.6) 
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first three columns of Table 1 are 5, 217 and 13108, and 3, 137 and 8150 for a centre gluon quark 
vertex. 



Green's function 


One loop 


Two loop 


Three loop 


Total 




3 


19 


282 


304 


c A c B 


1 


9 


124 


134 




1 


6 


79 


86 




2 


33 


697 


732 


Total 


7 


67 


1182 


1256 



Table 2. Number of Feynman diagrams for each Green's function for the Curci-Ferrari gauge 

renormalization. 



To deduce the renormalization constants themselves for each Green's functions, we apply 
the procedure discussed in 0. Here one computes the Green's functions in terms of the bare 
parameters, g , a Q , a Q and Co- The renormalized values are introduced by the definitions 
(J3.1|) and iteratively by loop order the renormalization constants are fixed by demanding that 
the overall infinity remaining is absorbed by the renormalization constant associated with that 
particular Green's function. As the Mincer algorithm is based on dimensional regularization 
in d = 4 — 2e dimensions, we have absorbed all the poles in e using the modified minimal 
subtraction scheme. 

From a practical computing point of view we organised the one and two loop renormalization 
in a different way from the three loop computation. For the former we retained an arbitrary a, 
a and ( in the extraction of the renormalization constants. However, for the three loop case, due 
to the increase in the number of actual algebraic terms in a Feynman diagram to be evaluated, 
due to the presence of the parameters in the propagators and vertices, we chose to fix ( to be 
or 1 and a = when the Feynman rules were substituted. This speeded up the computation 
significantly and avoided very large intermediate Form files which are generated. Running our 
code in the Landau gauge first allowed us to check the programme was performing correctly 
before generating the explicit value of the Feynman graph in the MAG. Moreover, at three loops 
one does not need to be concerned about the renormalization of the bare C i n this approach 
since any corrections to this would only appear at three loops. 

Having summarized the details of the computation we now record the explicit results. Rather 
than present the renormalization constants themselves, we have encoded them in the renormal- 
ization group functions defined by (|3.6j) . Hence, with a = £ = 0, we find, 

7a(«) = ^ [N% ((3a - 13)CU + STpNf) + N% ((-3a + 9)C A ) 
1 

+ 



1 

+ 



3456A2 



o3 



NX 2 ((6a 2 + 66a - 354) C\ + 2AQC A T F N f + l92C F T F N f 

+ N£N$ ((3a 2 + 210a + 331)C^ - 80C A T F Nf) 
+ Nf ((15a 2 - 6a - 33)Ci)l a 2 

A^ 3 ((162a 3 + 2727a 2 + 2592aC 3 + 18036a + 1944C 3 - 1 19580) C\ 

+ (- 6912a - 62208C 3 + Yimi2)C 2 A T F N s 
+ (82944C 3 + 960)C A C F T F N f - 29lMC A T F Nf 
- m\2C 2 F T F N f - 16896)C F T F N?) 



12 



+ NfN% ((2133a 3 + 162a 2 C 3 + 25785a 2 + 14904aC 3 + 61479a 

- 3564C 3 + 105550)Cl 
+ (- 13392a - 62208C 3 - Zl2M)C 2 A T F N f 
+ (82944C 3 - 77760)C A C F T F N f - 8960C A T F Nj?) 

+ iV|iv/((-324a 3 C 3 + 1728a 3 - 6480a 2 C 3 + 14256a 2 

- 11988aCs - 26298a - 129924C 3 - 113751)Cl 
+ (13392a + 41472C 3 + 9936)C A T F Nf) 

+ n/((- 4536a 3 Cs - 270a 3 - 18792a 2 Cs - 3294a 2 - 82296aCs 

+ 42714a - 176904C 3 + 101952)Ci)l a 3 + 0(a 4 ) (3.7) 

where Q n is the Riemann zeta function and 

7a(a) = j^TV? l N * ( (_ 3a ' + 26a)C<A " WaTFN f) + N A ((" 6 « 2 " 36a " 36 )Ca)] a 

NX 2 ((- 3a 3 - 51a 2 + 354a)Cl - 240aC A T F N f - 192aC F T F N f ^ 

+ NXNX 1 ((- 27a 3 - 339a 2 - 647a - 928)Ci + (160a + 512)C A T F N f ^j 
+ N f ((- 30a 3 - 366a 2 + 294a + 2016)C*l) 



12aNX 
1 

+ 



48aiVf 



+ 



6912aiVf 



a 2 



2 



iVf ((- 162a 4 - 3348a 3 - 5184a 2 C 3 - 25218a : 

- 3888aC 3 + 239160a)Ci 
+ (7344a 2 + 124416aC 3 - 349824a)ClT F 7V / 
+ (- 165888aC 3 - 1920a) C A C F T F N f + 58368aC A T F Nf 
+ 13824aC 2 ,T F iV / + 33792aC F T F Nf) 
+ NX 2 NX L ((- 2754a 4 - 48492a 3 - 14256a 2 C 3 - 155493a 2 

+ 27864aC 3 - 256744a + 209952C 3 - 548904) C\ 
+ (29376a 2 + 207360aCs + 36064a 
+ 331776C 3 + 136128)C A T F N f 
+ (- 331776aC 3 + 311040a 

- 663552C 3 + 705024)C A C F T F N f 
+ (35840a + 61440)CUT|iV / 2 ) 

+ NXNf{{- 7884a 4 - 133920a 3 + 76464a 2 C 3 - 151524a 2 

+ 517752aC 3 + 503388a + 1666656C 3 + 1014012)C^ 
+ (29376a 2 - 248832aC 3 + 5184a 
- 995328C 3 - 812160)C A T F N f ) 

+ iV~/((- 6480a 4 - 105840a 3 - 220320a 2 C 3 + 110700a 2 
- 784080aC 3 + 373032a 



- 1021248C 3 - 3148632)Q) a 6 + 0(a 4 ) . (3.8) 

For completeness we record the sum of the previous two anomalous dimensions partly to indicate 
the singular nature of this renormalization group function, but also because it corresponds to 
the renormalization of the gauge parameter itself from the convention we have used to define it. 
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We have 
7a(o) + 7a (a) 



C A 



Next, 



7 Ai (a) 



4aiVf 



a 2 N° A 



(4a 2 + 6a + 12) iVf 



+ 



C A 



48a7Vf 



NX 2 



+ 



6912aiV£ 3 



((3a 3 + 15a 2 )C A ) 

+ NXNi ((- 24a 3 - 129a 2 - 316a - 928)C A 

+ (80a + 512)T F JV/) 
+ n£ ((- 15a 3 - 372a 2 + 261a + 2016)C A ) 

iVf ((162a 4 + 2106a 3 + 10854a 2 )Ci - 6480a 2 C 2 A T F N f ) 



+ 7V| 2 iVf((1512a 4 + 324a 3 C 3 + 3078a 3 + 15552a 2 (3 
- 32535a 2 + 20736aC 3 - 45644a 
+ 209952(3 - 548904)Cl 
+ (2592a 2 + 82944aC 3 - 26464a 

+ 3317760? + 136128)CiT F iV / 
+ (- 165888aCs + 155520a 

- 663552(3 + 705024) C A C F T F N f 
+ (17920a + 61440)CUT|iV / 2 ) 

+ N£Nf((- 648a 4 C 3 - 4428a 4 - 12960a 3 Cs 

- 105408a 3 + 52488a 2 Cs- 204120a 2 
+ 15552a 2 Cs + 257904aC 3 + 275886a 
+ 1666656C 3 + 1014012)Cl 
+ (56160a 2 - 165888aC 3 + 25056a 
- 995328Cs - 812160)C A T F N f ) 

+ iv/((- 9072a 4 Cs - 7020a 4 - 37584a 3 Cs - 112428a 3 
- 384912a 2 Cs + 196128a 2 - 1137888aCs 



+ 576936a - 1021248Cs - 3148632)Cl) 



+ 0(a 4 ) 



^[4T F N f -llC A ]a 



(3.9) 



+ 



+ 



1 



1 

54 



- MC\ + 20C A T F N f + l2C F T F Nf \ a 2 

- 2857C1 + 2830C A T F N f + 1230C A C F T F N f 

- 2>lQC A T F Nf - l08C F T F N f - 2UC F T F Nf 
and we have checked explicitly that when a = 

7ai (o) = - j Ai (a) + 0(a 4 ) . 
For the ghosts and quarks we have 



+ 0(a 4 



(3.10) 
(3.11) 



7c (a) 



1 

Inx 



NX ((a - 3)C A ) + NX 1 ((-2a - 6)C A ] 
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+ 



96iV 



2 



iVf ((6a 2 - 6a - \m)C\ + SOCUTpiV/) 

JV£N^ ((- 42a 2 - 126a - M7)C 2 A + 160C A T F N f ^ 

j2 / / . - O — ,o \ 1 o 



+ 

+ 7v/ ((12a 2 - 588a + 510)C 

1 

+ 



6912iV? 3 



2 ) 



+ - 



iV£ 3 ((162a 3 + 1485a 2 - 2592aC 3 + 3672a - 1944C 3 - 63268)C^ 

+ (- 6048a + 62208C 3 + 6208)C^7>JV / 
+ (- 82944C 3 + 77760)C A C F T F N f + 89mC A T F Nf) 
+ NX 2 NX 1 ((1242a 3 + 10287a 2 + 8748aC 3 + 2565a 
+ 57996C 3 - 19184)Cl 
+ (- 5616a + 103680C 3 - 4:7632)C A T F N f 
+ (- 165888C 3 + 155520)C A C F T F N f + 17920CaT|A^ 2 ) 

+ NXNf{(- 1296a 3 Cs - 1836a 3 - 16200a 2 Cs - 68148a 2 

+ 140292aC 3 - 161730a + 617868C 3 - 258174)Cl 
+ (35424a - 124416Cs - 18U4)C A T F N f ) 
+ iv/((- 18144a 3 Cs + 864a 3 + 11664a 2 Cs - 89532a 2 - 191160aC 3 

+ 128304a - 21384Cs - 135972)Cl)l a 3 + 0(a 4 ) (3.12) 

7cl(a) = 4^? [NX((-a-3)C A ) + Ni((-2a-6)C A )]a 

NX 2 ((- 6a 2 - 66a - 190)C\ + 80C A T F N f ^j 

+ N%N£ ((- 54a 2 - 354a - 323)C A + lWC A T F Nf} 
+ Nf ((- 60a 2 - 372a + 510)C^)] a 2 

NX 3 ({- 162a 3 - 2727a 2 - 2592C 3 a - 18036a - 1944C 3 - 63268)Cl 

+ (6912a + 62208Cs + 6208)Cl T F N f 
+ (- 82944C 3 + 77760)C ' A C F T F N f + 89mC A T F Nf) 
+ NX 2 Ni((- 2754a 3 + 648C 3 a 2 - 28917a 2 - 4212C 3 a 

- 69309a + 37260C 3 - 64544) C\ 
+ (25488a + 103680C 3 - im72)C 2 A T F N f 
+ (- 165888C 3 + l55520)C A C F T F N f + l7920C A T F Nf) 

+ iV|iv/((- 7884a 3 + 22680C 3 a 2 - 84564a 2 + 97524C 3 a 

- 47142a + 433836C 3 - 56430) C\ 
+ (25056a - 124416C 3 - 18144) T F N f ) 

+ Nf((- 6480a 3 + 34992C 3 a 2 - 70092a 2 + 8424( 3 a 

+ 114912a + 77112C 3 - 161028)C^)1 a 3 + 0(a 4 ) (3.13) 



6912iV? 3 



and 



^ {a) = ~N F ~ a 
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+ 



AN F 



576N F T F N£ 



a 2 + 22a + 23)C A C F N F + (a 2 - 14a + 2)N£C A T F 



6C F N F - 8C F NfT F N F 



(684a 3 + (1296C 3 + 3528)a 2 + (22464C 3 - 14094)a 



+ 30240Cs - ^2U)C\C 2 F N 2 F 

+ (- 810a 3 + (-2376(3 - 1908)a 2 + (-40608C 3 + 42039)a 

- 63072C 3 + l0miQ)C 2 A C F N F N%T F 
+ (180a 3 + (IO8OC3 - 1080)a 2 + (18576C 3 - 23211)a 

+ 27864C 3 - 39132)C 2 1 A2 2 r 2 , 
+ (- 5472a - 24l28)C A C F N F N£N f T F 
+ (3024a + 5760)C A NfN f T F 
+ (6912C 3 - 20592)C A C F N F N£T F 
+ 1728C F N F N£N f T F + 12WC F N F N£Nfl% 

+ 86AC F N F N£T F ] a 3 + 0(a 4 ) . 



Finally, for completeness we record that the /3-function emerges as 



/3(a) 



f o A - \t f n, 



34 



C 



+ 



28mC 2 A T F N f - 2857C1 + 12?>QC A C F T F N f - 316C A T F Nj? 

i a 4 

- W8C F T F N f - 26iC F T F N^ — + 0(a 5 ) . 



(3.14) 



(3.15) 



To have confidence in the correctness of these results it is important to indicate the checks we 
have carried out. First, from a renormalization point of view, using the method of jS] the double 
and triple poles in e at three loops and the double pole at two loops in the renormalization 
constants are not independent of the previous order one loop poles. Therefore, we have checked 
that these emerge correctly for both the Landau gauge and the MAG. This is a non-trivial 
observation given the particular structure of the renormalization group functions which depend 
not only on a but also on the colour group Casimirs and for the MAG, the dimensions of the 
centre and off-diagonal sector of the Lie group. Second, the /3-function correctly emerges from 
the renormalization of the centre gluon. Again this is non-trivial since in the renormalization 
of the 2-point function its divergence has to emerge to be independent of not only a but also 
of and N£ as well as being equivalent to the actual /3-function itself. By the same token we 
can of course trivially record that the four loop anomalous dimension for A 1 is, [S], 



j A i{a) 



3 [4T F Nf 



11C A ] a 



+ 



+ 



+ 



54 



MC A + 20C A T F N f + Y2C F T F Nf 



2857C1 + 28WC 2 A T F Nf + 1230C A C F T F N f 



31QC A T F Nf - l08C F TpNf - 2MC F T F Nf 



44 



+ 



,<3 
/656 



150653 ^ 4 ^ 39143 
^86~ ) Ca + I ~8~T 



136 



C.3 



7073 \ 



243 J 



ClCpTpN f + 



3 

/4204 
V^27~ 



C 3 C A T F N f 



352 

— C3 ) C A CpT F N f - mC a F TpN f 



16 



/7930 224 \ t2 /704„ 1352\ - „ 



( 17152 448 



424 



+ ^-Cs C A C F T F Nf - ——C A T F Nf 



1232 



C F T F Nf 



+ 



V 243 9 / r 1 243 r J 243 

80 _ 704 \ dj BCD dj BCD /i 6 64 _ 512\ d F BCD dj BCD 
9 3 C 7 iV F + V 3 C3 9 ) N F f 



/704 512 \ d F 



^ABCD^ABCD 



N F 



-W 



a 5 + 0(a b ) 



where 



d F BCD 



= Tr (t a t( b T c T d ^ 
and <i^' BC * D is d F BCD evaluated in the adjoint representation. 



(3.16) 



(3.17) 



The next checks concern the anomalous dimensions themselves in certain limits. We have 
already indicated that the programmes we have used correctly reproduce all the Landau gauge 
results prior to switching to the MAG. However, the anomalous dimensions are also related to 
those of the Curci- Ferrari gauge. For instance, for the off-diagonal gluon, off-diagonal ghost and 
quark, taking the formal limit N^/N^ — ► 0, then the following anomalous dimensions arise for 
arbitrary a, 

lim lA {a) = \ [(3a - 13)C A + 8T F N f ] a 



+ 



+ 



48 



3456 



(6a 2 + 66a - 354) C\ + 24QC A T F N f + l92C F T F N f ~\ a 2 

(162a 3 + 2727a 2 + 2592aC 3 + 18036a + 1944C 3 - 119580)Ci 

+ (- 6912a - 62208(3 + 174912)Cl T F N f 
+ (82944C 3 + 960)C A C F T F N f - 29184C A T F Nj? 

- 69l2C F T F N f - 16896)C* F T|iV/l a 3 + 0(a 4 ) (3.18) 



lim 7 a (o) 
at£/jv°-o 



12 



[(-3a + 26)C A - 16T F N f ] a 



+ 



+ 



48 
1 



6912 L 



lim [7.4(a) + 7 a (a) 
nX/nx^o 



(- 3a 2 - 51a + 354)C\ - 240C A T F N f - 192C F T F N f j 

(- 162a 3 - 3348a 2 - 5184aCs - 25218a 
3888Cs + 239160)Cl 



+ (7344a + 124416Cs - 349824)C^7>./V f 

+ (- 165888C 3 - l920)C A C F T F N f + 58368C*AT|iV/ 

+ 13824CfT F iV / + 33792C F T|Ar/l a 3 + 0(a 4 ) (3.19) 



aC A a(a + f>)C 2 A 9 

a H —a 

4 16 

3a 



128 



a 2 + 13a + 67)C A - 40T F iV / ] C A a 3 + 0(a 4 ) (3.20) 



lim 7c (a) 

nX/nx^o 



-[(a-3)C A ]a 



17 



+ 



1 

96 
1 



6912 



(6a 2 - 6a - 190)C| + 80C A T F N f 

162a 3 + 1485a 2 - 2592aC 3 + 3672a - 1944C 3 - 63268)C^ 
+ (- 6048a + 62208C 3 + 6208)C%T F N f 
+ (- 82944C 3 + 77760)C A C F T F N f + 8960C A T F N? 



and 



lim 7,/, (a) 



+ 0(a 4 ) 



aC F 1 r 
4 4 



(3.21) 



3a + 25)C A C F - 6C| - 8C F T F N f 



+ 



1 

288 



(27a 3 + 270a 2 + 216aC 3 + 2367a - 2484C 3 + 18310)ClC F 

+ (- 1224a - 9l84)C A C F T F N f + A32C F + 8UC F T F N f 

+ (3456Cs - W296)C A C F + 640C F T F N?)\ a 3 + 0(a 4 ) (3.22) 



where to take the limit for the quark anomalous dimension* we have used the result that 



n 1/ n x^° Np 



C F . 



(3.23) 



Comparing these limits with |121 1511 I52j . we observe that they are equivalent to the three 
loop MS anomalous dimensions in the Curci- Ferrari gauge for arbitrary a. That this result 
appears is not unexpected since Kondo indicated in 19. that the off-diagonal sector is in fact 
the Curci-Ferrari gauge. Indeed this observation, and its relation to the generation of a non- 
zero vacuum expectation value for the operator O, was one of the reasons for the recent renewed 
interest in both the Curci-Ferrari gauge and MAG. That the Curci-Ferrari anomalous dimensions 
correctly emerge is an important check on the full MAG computation. A final more trivial check 
rests in taking the formal abelian limit in the Landau gauge, C A — ► 0, C F — ► 1, T F — > 1 and 
a — * 0. One observes that both ghost anomalous dimensions vanish, the centre and off-diagonal 
gluon anomalous dimensions reduce to the quantum electrodynamics /3-function and the quark 
anomalous dimension tends to the electron anomalous dimension. 

Having justified the results for the full renormalization of the QCD Lagrangian in the MAG, 
we can now deduce the anomalous dimension of O. At two loops we actually computed the 
anomalous dimension directly by the same method as The operator was inserted in an 

off-diagonal ghost 2-point function and the corresponding renormalization constant Zq was 
extracted. Computing the associated anomalous dimension directly from the renormalization 
constant, the resulting two loop value correctly satisfied (|3.3|) . At three loops we took the point 
of view that the three loop 7 c i (a) was correctly determined and therefore used 1)3. 3 j) to deduce 



70 (a) 



1 



12N£ 



+ 



96N£ 2 



NX ((- 3a + 35)C A - 16T f N f ) + N% ((- 6a - 18)C A 

NX 2 ((- 6a 2 - 66a + 898)C\ - 560C A T f N f - 384C F T f N f 

+ N£N% ((- 54a 2 - 354a - 323)C^ + l60C A T f N f 
((- 60a 2 - 372a + 510)0^ 



+ Nf 



*Whilst the renormalization constant for Z$ was explicity given in |12| for the Curci-Ferrari gauge, the actual 
anomalous dimension was inadvertently omitted. 
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+ 



6912iVf 



NX ((- 162a d - 2727a 2 - 2592C 3 a - 18036a - 1944C 3 + 302428)Q 

+ (6912a + 62208C 3 - 356032)C\T F N f 

+ (- 82944C 3 - 7mm)C A C F T F N f + AM08C A T F Nj? 



13824C F T F N f + 337§2C F T F Nf) 



+ NX N%((- 2754a 3 + 648a^3 - 28917a 2 - 4212aC 3 
- 69309a + 37260C 3 - 64544) C\ 
+ (25488a + 103680C 3 - 13072)CiT F iV / 
+ (- 165888C 3 + 155520)CU C F T F N f + 17920C A T F Nj?) 

+ N£Nf((- 7884a 3 + 22680a 2 C 3 - 84564a 2 + 97524aC 3 - 47142a 
+ 433836C 3 - 56430) C\ 
+ (25056a - 124416C 3 - 18144)C^Z>JVf) 

+ Nf((- 6480a 3 + 34992a 2 C 3 - 70092a 2 + 8424a<C 3 + 114912a 



+ 77112C 3 - 161028)C|) a 3 + 0(a 



(3.24) 



Again there are several checks on this result aside from the internal renormalization group 
consistency check. First, we verified that the Landau gauge anomalous dimension emerged 
correctly. Second, in the formal limit N^/N^ — > 70(a) tends to the Curci- Ferrari gauge result 
for all a, H2EDE2]. 

Finally, having derived the anomalous dimensions of all the fields and O for an arbitrary 
colour group, we record the explicit results for the two main Lie groups of interest. For SU(2) 



with C A = 2,C f = 3/4, T F = 1/2, N F = 2, Ni 



2 and N% 



1, we have 



7a(o) 



7« (a) 



1a{o) + 7« (a) 



[3a - 17 + 4N f ] - + [45a 2 + 678a - 787 + 272N f 



+ 



48 

2592a 3 C 3 + 6507a 3 - 15552a 2 C 3 + 75087a 2 - 10260aiV) 
12960aC 3 + 190161a - lOOOOiV/ - 36288iV/C 3 + 146572iV) 



+ 
+ 



+ 
+ 



- 217728C 3 - 329995] 
3a 2 + 4a - 9 - 2aN f 



1728 

a 



+ 0(a 4 ) 



3a 



a 

6a 



- 12a 3 - 156a 2 + 52a + 20 - (29a - 32) Nf 

- 2160a 4 - 37152a 3 + 3672a 2 iV) - 10368a 2 C 3 - 63504a 2 

+ 2780aiV/ + 5184aiY/C 3 - 33869aiV) + 20736aC 3 + 141632a 



+ 960JV/ - 
a 2 + 3a + 6 



25920iV/C 3 - 348iV/ + 196992C 3 - 207264 
a 

2a" 



432a 



+ 0(a 4 



51cT - 570a 2 



a 



48a 



371a + 160 + (40a + 256)^/ 

- 2592a 4 C 3 - 2133a 4 - 15552a 3 C 3 - 73521a 3 + 4428a 2 A?/ 

- 54432a 2 C 3 - 63855a 2 + mOaiV/ - 15552aA^/C 3 + 11096aA^ 

- 134784aC 3 + 236533a + 3840JV/ - 103680JV/C3 - 1392./V/ 
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+ 787968(3 - 829056] 



1728a 



+ 0(a 4 ) 



lAi («) 



7c (a) 



7c* (a) 



7o(a) 



2[A^ - 11] | + [497V/ - 272]^- 



+ 



16607V/ 2 + 5241 77V,- - 182848 



432 + °<°> 



3a a + 



3a 2 - 54a - 59 + 107V/ 



l a 



+ - 648a 3 C 3 + 108a 3 - 648a 2 C 3 - 5400a 2 + 3240aC 3 - 4860a + 2807V/ 2 
+ 12967V/C 3 + 22617V/ + 44712C 3 - 38600] — + 0(a 4 ) 

. „2 

- [a + 3] a + - 6a 2 - 42a - 28 + 57V/ 



+ 



1080a 3 + 2592a 2 Cs - 11772a 2 + 1620a7V/ + 5184aCs - 12528a 

i a 3 

+ 2807V/ + 12967V/C 3 + 33417V/ + 33696C 3 - 32444 — + 0(a 4 ) 



7v(a) = -a a + 



4a 2 + 152a + 265 - 24iV/ 

i c;qqq„.2 



l a 
32 



+ 0(a 3 ) 



672a' 5 + 576a^C3 + 5328a^ - 1728a7V/ + 10944aC 3 + 10848a 

i a 3 

+ 160/V/ - 98207V/ + 6912C 3 + 50863 — + (9(a 4 ) 



[- 3a + 13 - 27V)] - + 



384 

4a 2 - 28a + 72 - 137V, 



a 

y 



+ 



720a 3 + 1728a 2 C 3 - 7848a 2 + 1080aTV/ + 3456aC 3 - 8352a 



+ 7407V/ + 8647V/C 3 - 152457V/ + 22464C 3 + 39320J — + 0(a 4 ) . 

(3.25) 



The one loop expressions agree with the limited known results. Repeating the same exercise for 
SU(3) with C A = 3,C f = 4/3, T F = 1/2, 7V> = 3, 7V£ = 6 and TVj* = 2, we have 



7A(a) = [3a - 15 + 27V/] - + 39a 2 + 609a - 1113 + 2247V/ 



7a (a) 



l a 
24 



- 1836a 3 C3 + 9495a 3 - 12258a 2 C3 + 115056a 2 - 14832a7V/ 
+ 28620aC 3 + 334701a - 99207V/ 2 - 437767V/C 3 + 2297047V/ 



- 182088C3 - 839337] 



1152 



+ 0(a 4 ) 



15a 2 + 42a - 36 - 8a7V/ 



12a 



+ 
+ 



48a 



- 138a 3 - 1842a 2 + 1539a - 768 - (408a - 256)7V/ 

- 59292a 4 - 1034424a 3 + 91800a 2 7V/ - 259200a 2 Cs - 2424195a 2 
+ 64000a7V/ 2 + 193536a7V/C 3 - 1301712a7V/ + 914976aC 3 

+ 6029820a + 153607V/ 2 - 4423687V/ ( 3 + 2681287V/ 



+ 5868288C 3 - 5046732] 



a 

6912a 



+ 0(a 4 
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1a{o) + 7a (a) 



7a* (a) 



7c (a) 



7 c <(a) 



a 2 + 6a + 12 



a 

4a 



+ 
+ 



- 60a 3 - 624a 2 - 687a - 768 + (40a + 256)7V/ 



48a 



+ 0(a 3 ) 



11016a 4 Cs - 2322a 4 - 73548a 3 Cs - 344088a 3 + 2808a 2 7V/ 

- 87480a 2 C 3 - 415989a 2 + 4480a7V/ - 69120a7V/C 3 + 76512a7V/ 

- 177552aC 3 + 993798a + 153607V/ - 4423687V/C 3 + 2681287V/ 



+ 5868288C 3 - 5046732] 



a 



6912a 



+ 0(a 4 ) 



W-33]| + 2 [197V/- 153] y 



+ 



- 3257V/ + 150997V/ - 77139] ^ + 0(a 4 ) 



[a - 15] - A + 



a 
96 



+ 



60a^ - 1020a - 2241 + 2007V/ 

- 22032a 3 C 3 + 10908a 3 - 36936a 2 C 3 - 161298a 2 - 17928a7V/ 
+ 238464aC 3 - 234657a + 112007V/ + 967687V/C 3 + 2066167V/ 



+ 2301696C 3 - 2791386] 



6912 



0(a 4 ) 



- 5 [a + 3] | + [- 276a 2 - 2028a - 2169 + 2007V/] ^ 

+ - 59292a 3 + 108864a 2 Cs - 657666a 2 + 81864a7V/ + 193104aCs 
- 1137267a + 112007V/ + 967687V/ Cs + 2584567V/ 



+ 1661472(3 - 2619450] 



6912 



+ 0(a 4 



7v( a ) = 



a a + 
+ 



- 3a 2 + 138a + 262 - 167V/ — + 0(a 3 ) 



8532a 3 + 5832a 2 C 3 + 71496a 2 - 19224a7V/ + 117936aC 3 + 210195a 

i a 3 

+ 12807V/ - 1142407V/ +43848C.3 + 948012 — — + 0(a 4 ) 



7o (o) = [- 15a + 87 - 87V/] 



12 



+ 
+ 



a 
96 



- 276a^ - 2028a + 7623 - 10167V/ 

- 19764a 3 + 36288a 2 Cs - 219222a 2 + 27288a7V/ + 64368aCs 

- 379089a + 176007V/ + 322567V/Cs - 5580727V/ 



+ 553824C 3 + 2418114] 



2304 



+ 0(a 4 



(3.26) 



4 Discussion. 



We have provided a comprehensive discussion on the three loop MS renormalization of QCD 
in the maximal abelian gauge. Indeed this article represents the first calculations beyond one 
loop as well as the first for Lie groups other than just SU(2). By explicit computation we 
have determined all the anomalous dimensions and /3-function before deducing the anomalous 
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dimension of O at three loops in MS. Indeed it is the explicit expression for the latter which 
will be the key to studies of the condensation of O in the MAG which we hope to examine next. 
One useful observation from the main results is the relation of the MAG anomalous dimensions 
to those of other gauges and in particular the Curci-Ferrari gauge. That the results for the 
latter appear in the formal limit N$/N£ -> is reassuring, though their prior existence was also 
of a more practical use in helping to establish the veracity of the final MAG renormalization 
group functions. Though from the actual structure of the final expressions it is clear that they 
could not be constructed from knowledge of the same anomalous dimensions in the Landau or 
Curci-Ferrari gauges. 
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A Feynman rules. 



In this appendix we record the Feynman rules we used for the maximal abelian gauge fixing 
in momentum space which are derived from (|2,5|) and (|2.13|) using a symbolic manipulation 
programme written in FORM. For the propagators we have 





-p)) = 




-p)) = 


(c a (p)c b (- 


-p)) = 


<cW(- 


-p)) = 




-p)) = 



5 ab 

ft 
p 2 



Vpu 
Vpu 



(1 



ft 



(1-a) 



PpPu 
PpPv 



P 



6 ab 
p2 
§ij 
p2 

P 2 



where p is the momentum. The non-zero 3- and 4-point vertices are 

(^(P1$(P2)^(P3)) 

(4>i)^2)<Kp 3 )> 
(A a ( Pl )c b ( P2 y(p 3 )) 



i 



TP! - P3 



(A*(pi)c b (p 2 )c k (p 3 )) 
(A«( Pl )cJ( P2 )c c (p 3 )) 
(A i ll (p 1 )c b ( P2 )c c (p 3 )) 
{Al{ Pl )A h AP2)Al{p 3 )) 
(Al( Pl )A b u (p 2 )Al(p 3 )A d p (p 4 )) 



(A.l) 



-igf abc K 9 

- igf abk (-Cps), 
-igf acj (pi+P3), 
-igf bci (-Pi-2p 3 + P3 C), 

ig jabc ( Vua ( p2 _ + ^(pg - pi ) v + ^(pi - p 2 ) a ) 

ff Cd {-V^Vup + VwVva) + fT M {-VvvVap + Vl*pVv<r) 



adbc 
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<A«(pi)4(P2)A£(p3)4(P4)> 

(A^ Pl )A b u (p 2 )A k MA l p ( Pi )) 



(A;( Pl )A b u(P2)c c (p3)c d (p 4 )) 

(A^ Pl )Ai(p 2 )cf(p 3 )c d (pt)) 

(A^ Pl )Al(P2)c c (p 3 )c l (p 4 )} 
(A i fl {p 1 )Al(p 2 )cf{p 3 )c d (p,)) 

(c a ( P l)c b (p 2 )^(p3)c d (p 4 )) 
(?(pi)c b (p2)c C (p 3 )c l (p 4 )) 

(c a ( Pl y(p 2 )c c (p 3 )c l (p 4 )} 

where the momentum flow for each field is into the vertex. We note that we have recorded the 
Feynman rules as generated from the full MAG Lagrangian, using a Form routine, without 
recourse to the simplifying properties of the Jacobi identity of the Lie algebra. For example, 
the final rule of the non-zero set, QA.2JI . actually vanishes after application of a Jacobi identity. 
However, in the construction of the routines to perform the overall calculation, we have relegated 
all the algebra associated with the group theory to a common Form module which encodes the 
necessary simplifying lemmas. That module is placed after the module where the above Feynman 
rules are substituted. The remaining Feynman rules for the 3- and 4-point vertices which have 
not been recorded above are trivially zero due to the fact that they would involve either two or 
more centre indices in the case of the 3-point vertices or three or more centre indices for the 
4-point vertices. 
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